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Abstract 

We consider a laser induced molecular excitation process as a decay of a sin- 
O gle energy state into a continuum. The analytic results based on Weisskopf- 

^3 



Wigner approach and perturbation calculations are compared with numerical 
wave packet results. We find that the decay model describes the excitation 



Qh! process well within the expected parameter region. 

1 Introduction 

It has been known for a long time that when an isolated bound state is interacting with a 
continuum of quantum states, the occupation of the bound state experiences exponential 
decay; the first proper treatment was given by Weisskopf and Wigner in 1930 p]. A simple 
example of such a situation is a two-state atom coupled to the electromagnetic modes of 
the vacuum, which leads to spontaneous emission; see e.g. refs. [Q, The general model 
has been described in detail within the formalism of scattering theory in ref . [f| , but this 
approach provides only a steady state treatment. In a recent paper || we considered this 
problem in a genuine time dependent setting and discussed the decay process as a method 
to prepare a moving wave packet on a molecular electronic state in accordance with our 
earlier treatment of coupled channel molecular wave packet dynamics; for a review see 
ref. |J. Now we proceed to study the validity of this exponential decay description of 
laser induced wave packet excitation in molecules. 

The prototype model for the Weisskopf-Wigner decay is shown in fig. |T](a); the discrete 
state embedded in a continuum acquires a Lorentzian line shape, which corresponds to 
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exponential decay of its occupation. The theoretical description leading to exponential 
decay is based on an approximation which represents the continuum spectral density by a 
single dominating pole 0, |J . In the general case this is only an approximation, which is 
expected to hold in the perturbative limit of weak coupling between the discrete state and 
the continuum, in which case the result of time dependent perturbation theory agrees with 
the Weisskopf-Wigner result. When we go beyond this limit, the result is nonexponential 
decay M or a continuous transition of the decay into a periodic oscillational behaviour of 
the Rabi type as discussed in ref. ||. Thus we need to establish for molecular excitation 
the parameter regions where a) the perturbation approach is applicable, and b) where the 
excitation process can be described as exponential decay into the continuum. 




Figure 1: (a) An isolated energy level at E is coupled to a system with a continuous energy 
spectrum E. (b) In the molecular excitation model the isolated state is the ground state 
of the molecule, Ui, and the continuum is provided by the linear excited state potential 
U 2 . 

In this paper we develop further the decay model introduced in ref. M and featured 
in fig. |l](b). The lowest vibrational level of the molecular electronic ground state corre- 
sponds to the discrete state; it can be described adequately by the harmonic potential U\. 
The molecular counterpart of the continuum is an excited electronic state which can be 
adequately described by the linear potential U% near the molecular equilibrium position. 
Usually state 2 is a dissociating one, but it could also represent a bound state in a energy 
region where the vibrational states have a very small separation (quasicontinuum). These 
potentials are a good basis for the molecular processes within the Born-Oppenheimer 
approximation. If we also apply the dipole and rotating wave approximations to the 
laser-molecule interaction, we can shift the excited state potential down by one laser pho- 
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ton, which then leads to the elimination of the rapidly oscillating field term from the 
coupling between the states 1 and 2; for more details see ref. |J. 

If the parameters of our model are chosen such that the vibrational excited states of 
the harmonic potential remain mostly unoccupied during the interaction, the decay of 
the oscillator ground state can be taken to be a model of the level scheme in fig. [I]. The 
couplings between the oscillator state and the eigenstates of the slope are not constant; 
the Condon factors describing the overlap between the ground and excited state wave 
functions have an energy dependence which complicates the model. Thus for simplicity 
(one parameter less to consider) we arrange the sloping potential Ui so that a Franck- 
Condon transition at the origin of the coordinate system gives the maximum overlap 
between the wave functions. Then we expect the Condon factor to depend only slowly 
on energy around E and the model will simulate the behaviour of the Weisskopf-Wigner 
situation. 

This paper is organised as follows. Section ^| presents the basic model. In sec. |3| we 
give the results of our numerical computations and discuss the occurrence of exponential 
decay. It is found that this can be observed in a parameter range where the decay rate 
is less than a definite value. This indicates that the exponential behaviour can indeed be 
achieved only in a perturbative regime as expected. We also find that the occupation of 
the state 2 emerges in the form of a localized wave packet even beyond this perturbative 
parameter regime. In sec. [| we discuss the analytic result for the decay and compare this 
with our results obtained from numerical integration of our wave packet model. Finally 
we summarize and discuss our results in sec. |5|. 



We consider a one dimensional system, spatial coordinate x, with two internal states 
{|1), 1 2)}. The Hamiltonian of the system is written as 



where the coupling V is chosen real. We define (in suitably chosen dimensionless units, 



2 The model 



H=[T X + U x (x)\ |1)(1| + [T x + U 2 (x)} |2)(2| + V [|2)(1| + |1)(2|] , 



(1) 



cf. refs. § 1) 



T 



X 



dx 





(2) 
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The Hamiltonian ([!]) describes our chosen model as shown in Fig. |T](b). Since in our 
scaling the energy of the lowest level of the harmonic potential is 

Eo = r = 7i' (3) 

the energy of the linear potential resonates with the ground state at x = 0, see fig. [p. 
The eigenfunction of the ground state of the harmonic oscillator is 

V?o(*) = ^^exp(-* 2 /2 3 / 2 ) (4) 

with the energy eigenvalue (|3|). In the domain we are investigating, the higher eigenfunc- 
tions are not supposed to play any role, and they are consequently not needed, as already 
mentioned in sec. [I]. 

The eigenvalues of the linear state Ui form a continuum {—00,00} with the corre- 
sponding eigenfunctions easily obtained by a Fourier transform 

Mk) = Jdxe- ikx ij E (x). (5) 

The functional form is that of the Airy function, and a suitable normalization is given by 

(6) 



^Pe(x) = I dk exp 



2k ^[a 



1 [ x H — fc — 1 — 



otl a 



where the energy parameter is 

e-E-j,. (7) 

It is easy to prove the energy normalization 

J %l)* E ,(x)tl) E (x)dx = 5{E-E'). (8) 

and the completeness 

r ij* E (x')ijE(x) dE = d(x-x'). (9) 



With these relations we see that the density of states is unity. 
We assume that the system is prepared in the initial state 

|z>=y?o(x)|l>, (10) 

which is then coupled to the continuum by the parameter V. Using the result of the time 
dependent perturbation theory (Fermi Golden rule) the leakage into state |2), i.e. into 
our final state |/), is expected to occur at the rate 

r = 27r|(z|iJ|/)| 2 = 27ry 2 |^ |^ )| 2 , (11) 
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where the continuum function is evaluated at the energy (|3|) and the parameter dependence 
is almost entirely in the Condon factor (^olV'-Eo)- 

From the normalization of eq. (|) we expect the decay rate to scale as 

V 2 

r oc 2tt — = r , 12 

a 

where we have defined an effective decay parameter r . When T is below some limiting 
value, we expect that eq. (|ll|) gives a reasonable approximation to the decay rate of the 
initial state (|T0|). The result is modified by the Condon factor which will be discussed 
below. 



3 Numerical calculations 

We have used the numerical approach described in sec. 2.3.3 of ref. || to propagate the 
initial state wave packet (|T0|) on the coupled energy surfaces of eq. (p]). However, instead 
of switching the coupling on suddenly we have used the coupling term 

[ V sech[(t-t )/T)], t<t 
V(t) = (13) 
[V t>t 

where t indicates in practice the beginning of the excitation process, and we have set 
T = 0.05 in our scaled units. 

We have especially looked into the possible occurrence of exponential decay and the 
emergence of an outgoing wave packet on state |2) in our model. Some typical results of 
our calculations are shown in figs. ||] and |3|. They present the logarithm of the occupation 
Pi of the state |1) as a function of time. Both figures show descending straight lines 
which indicate exponential decay. For very large values of r the initial state occupation 
Pi becomes clearly oscillatory as demonstrated in fig. |j. 

There is no absolute criterion when exponential decay is observed; the occurrence of 
linear behaviour in the semilogarithmic plots is subject to an arbitrary criterion. Running 
various parameter combinations, however, indicates that monotonic exponential decay 
is observed for values of V 2 smaller than a constant times a, thus indicating that the 
parameter plays a decisive role in the phenomenon. In fig. ^| we show the data 
points that correspond to exponential and nonexponential decay in the (a,V 2 ) plane. 
The exponential decay points fall into a region for which we have roughly 

r <2...2.5, (14) 
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t-tr 



Figure 2: The occupation P\(t) for Vq = 0.5. Note the logarithmic scale on the vertical 



axis. 



which is thus taken to be the limit of validity of the simple pole approximation for decay. 
Typically the nonexponential behaviour appears in P\ and P2 as oscillations, and the 
exponential behaviour as a steady change; this is demonstrated in fig. |[ 

We have also studied how the emerging probability on level |2) forms a wave packet. 
For the case V = 0.7, a = 2 (r — 1.5) we observe a clear wave crest emerging in fig. [7|. 
This is followed by a long, slowly diminishing tail representing the final leakage out of state 
We have, however, observed that the formation of a wave packet is not conditioned on 
the occurrence of exponential decay. In fig. |8| we still see a well developed wave packet for 
the parameters Vq = 1.5 and a = 3 (r ~ 4.7) far exceeding that for which perturbation 
theory would be expected to work. 

We have not carried out a systematic investigation of the occurrence of a well formed 
wave packet. This is as subjective a phenomenon as the exponential decay. Because we 
initially find no probability on state |2) and ultimately there is no longer any probability 
leaking out of state |1), there must be a localized wave packet emerging out of the system. 
However, a few words connecting the results of this paper to our previous work pi fill 



l!g, [13[] are appropriate here. 

In both figures |7| and |8| we see a clear, peaked wave crest wave packet. On the other 
hand, figures ^| and |3| show that for all parameter values the ground state occupation 
behaviour is not exponential when t ~ to- This is because initially the excitation process 
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follows the well known area theorem, which for two resonant levels and a coupling Vq 
states that P x = 1 - P 2 = cos 2 [K (t - to)] - 1 - Va(t - t ) 2 + 0[(t - t ) 4 ] for t > t . 
This is the origin of Rabi oscillations (if the levels are off resonance, then the oscillation 
frequency and amplitude are modified). For more details see e.g. ||. These oscillations 
simply indicate that both states are "equal", i.e., transitions can take place also back 
from state |2) to state |1). In the Weisskopf-Wigner regime the second state acts like a 
reservoir and population can move only from state |1) to state |2). 

This is the point where the dynamics of the wave packet on state |2) comes to the 
rescue. In our earlier studies we have given special attention to cases where the dynamics 
inhibits the Rabi oscillations, because the excited state population is accelerated away 
from the resonance region faster than the Rabi oscillations take place [ |TU| , |TT| , |T2"| , |T3 |. 
Although our earlier studies were made for short pulses, the basic principle is the same. 
Thus the reservoir model works best for the excitation when a is large and wave packet 
motion on state \2) takes place swiftly after excitation. 

We have now established a loose criterion for obtaining exponential decay. Next we 
need to consider the validity of the perturbation result flTTD . The exponential behaviour 
is only a necessary but not a sufficient condition for eq. ([□]) to hold. For suitably large 
values of Vq the perturbation result may fail even if the criterion (|TJ) is fulfilled and the 
decay is exponential. For this purpose we calculate the full perturbation result for our 
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Figure 4: The occupation Pi(t) for Vo = 3 and a = 2. Here T — 28 and the oscillatory 
behaviour is obvious. 

model in the following section. 



4 Analytic considerations 

The perturbative rate ([11]) contains the Condon factor 

Soi = (Vo\lpE ), 



(15) 



which is the overlap of the initial and final states; see e.g. refs. |L4], [15], [L6[ |17| . This 
scalar product is most efficiently calculated in momentum space, where we need the 
Fourier transform of the state (|4]). Calculating the continuum wave function ipE for the 
energy (Q) we set £ = 0, which gives 

25/8^1/4 



5, 
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where we have used the substitution k = 

Ai(x) = 



du exp 



era u J 
~2 + 3a, 



la 



j\[2 — u. Using the Airy function 



ds 
2^ 
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Figure 5: The regions of exponential (o) and nonexponential (x) decay in the («,Vq 2 
plane. The shaded region corresponds to the upper limit given in eq. flU]). 

we can write (after setting u = a 1//3 s in eq. (|16D) 



25/8^1/4^/3^2 

<S>oi = — Al 



a 



1/6 



The result ( |i8|) seems to contradict the simple dependence on a from (0), but the asymp- 
totic relation Ai(x) ~ x _1 / 4 exp(— 2x 3//2 /3)/a/47t, valid for large x, restores the correct 
scaling with a in the perturbative limit. In fact, the asymptotic result for the Franck- 
Condon factor becomes 

lim \Sm I 2 



This is exactly what the Condon reflection principle would predict [[y], 16, 18]. According 
to it 



lim ISVi '"' 



'oil 



d(U 2 - U x ) 



dx 



\(p (x = X )\' 



(20) 



X=XQ 



where x is the point where ^(^o) = ^lO^o); i n our case x o approaches zero as a — > oo, 
and our ground state wave function contribution becomes \(fo(x = 0)| 2 = {2-n 2 )~ 1 ^. The 
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Figure 6: The time evolution of the excited state occupation P2 for a = 2 for various 
values of the coupling Vq. The transition from the steady behaviour to the oscillatory one 
with increasing Vq is clear. Here t = 1.0. 



reflection principle simply means that the excited state wave functions for steep potentials 
are so localised to x = x that they behave like S functions inside the overlap integrals. 
Thus we get that 

1 



lim T 

a— >oo 



(27T 2 ) 1 /' 



(21) 



The interesting conclusion from eq. flTED is that the decay rate is not the simple function 
of a suggested by eq. ( |T2] ) except in the asymptotic limit. The behaviour of however, 
is found to be simply monotonic with increasing a as shown in the fig. ||. 
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Figure 7: The excited state wave packet P2(x,t) for V = 0.7 and a = 2. This is in the 
exponential decay regime since To = 1.5. Here to — 1-0. 



5 Discussion 



We have shown that in the perturbation limit and for steep excited state potentials the 
Weisskopf-Wigner model can describe adequately laser induced molecular excitation pro- 
cesses. Our numerical results are in good agreement with analytical perturbation theory 
results in the parameter region where we expect them to be applicable. 

Although we have considered mainly the molecular excitation process, there are some 
interesting atomic systems that resemble the molecular model shown in fig. 0(b). One 
can trap single ions electromagnetically, and cool them down to the lowest motional state 
of the trap potential, which in this energy range corresponds to a harmonic potential |19 



ED , |2"T| , |2"2|j . Furthermore, one can also trap neutral atoms with magnetic fields, and cool 
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0.5 -i 




Figure 8: The excited state wave packet P2(x,t) for Vq = 1.5 and a = 3. This is in the 
nonexponential decay regime since Tq = 4.7. 



them evaporatively into the density-temperature region where Bose-Einstein condensation 
takes place |23|, [24j. The behaviour of the condensed atoms can be described as a single 
wave function 



25, 26 



Our model, therefore, resembles the recent studies for output 
couplers for the condensates, in which rf fields couple the trapping state to a nontrapping 
(continuum) state — such a system can be regarded as a first generation atom laser |27|, ^8 



Thus our model may provide some insights when one wishes to move from pulsed atom 
lasers towards cw atom lasers. 
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Figure 9: The a dependence of the decay rate V. We have plotted the wave packet 
results (open circles), the full perturbation result (solid line) [eqs. ([TTJ) and flT8|)] and 
the asymptotic result (dashed line) [eq. (pip]. All selected data points correspond to 
exponential decay 
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